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Approximate formulas for the inband RCS of an array with series feed 
have been (terived The formulas are based on the hypothesis that an incident 
wave excites forward and backward traveling waves on the main line. The 
approximate fmmulas are in good agreement with results obtained using 
scattermg matrices, thereby verifying the ass' /^ons made in the 
a;^oximate solution Spikes in the RCS pattern ikiv' i een identified with 
specific scattering sources in the array. The parameters affecting the level and 
location of the lobes have been noted The main advantage of tLe approximate 
method relative to the scattering matrix method is its speed. The scattering 
matnx method requires that a matrix equation be solved, and its size increases 
with the number of array elements. 
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I. INTRODUCTION 


Current state-of-the-art radars are sensitive enough to track small arms 
fire and even a swarm of insects. In response to this sensitivity, a higher 
standard of stealthy design has been employed to avoid detection by radar. 
The result has been radically stealthy platforms like the F-117A fighter and 
B-2 bomber. These are designed to have the smallest possible radar cross 
section (RCS), low infi'ared (IR) emissions, and low microwave emissions in 
order to avoid being detected and tracked. 

The RCS of a target not only depends on the physical shape and 
material, but also on its sub-components such as antennas and other sensors. 
These components on the platforms must also be designed to meet low RCS 
requirements in addition to their sensor system requirements. In some cases 
the onboard sensors can be the predominant factor in determining a 
platform’s total RCS. A typical example is a high gain antenna on a low RCS 
platform. If the antenna beam is pointed toward the threat radar and the threat 
fi-equency is in the anteima operating band, antenna scattering can be 
significant. 

Scattering fi’om antennas has been the subject of interest since the 
1950s. Extensive work has been done with regard to the determination of the 
antenna parameters [Refl], and dipole scattering and the effect of the 
terminal load impedance [Ref2]. The RCS of boms [Ref3], reflector 
anteimas [Ref 4], and microstrip elements and arrays [Ref S] have also been 
extensively studied. Detailed analysis of the inband scattering characteristics 
of arrays with parallel feed networks have also been investigated [Ref 6]. 
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This thesis examines the scattering properties of phased arrays with 
series feed networks (Figure 1) which operate at the same frequency as that 
of the illuminating threat radar. For this so-called inband case, the threat 
signal will penetrate into the feed and be reflected at the mismatches and 
junctions of the array. This occurs even for a well-matched array because 
small mismatches exist within the feed The total scattered power is a small 
fraction of the incident power, but for a large array, a large number of 
reflections exist which significantly contribute to the total RCS of the array. 
High performance arrays must not only meet the system operating 
requirements (i.e., high gain, low sidelobe levels, etc.) but also the RCS 
requirements. To optimize the tradeoff between antenna RCS and radiation 
performance, an efficient and accurate model of RCS is crucial. This thesis 
presents an approximate inband scattering mode! for arrays with series feeds. 
The approximate model is shown to be in good agreement with a more 
rigorous scattering matrix approach. The behavior of the RCS for various 
feed parameters is also determined. 

Chapter 11 presents the theoretical background for the RCS analysis of 
the series-fed array antenna, and contains the derivation of the formulas for 
the approximate solution. Chapter III presents some computed data for both 
the approximate and rigorous solutions. Results me aho presented to illustrate 
how the RCS is influenced by various feed parameters. Finally, Chapter IV 
summarizes the results of the data acquired and concludes with some design 
guidelines for controlling the behavior of the RCS. Some conclusions and 
recommendations are presented. 
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Pigarc 1: Series feed antenna network. 




II. A;4TENNA AND FEED SCATTERING ANALYSIS 


Scattered fields from a body are set up by a combination of induced 
electric and magnetic currents. The cuirents are induced such that the fields 
satisfy the boundary conditions and Maxwell’s equations. The scattered field 
can be determined from the current using the radiation integrals. However, in 
most practical cases, determining die induced current is difficult or 
impossible. For a target such as an antenna the prediction of RCS is difficult 
because of many scattering sources both at the aperture and in the feed 
circuit. In fact, even antennas with identical amplitude and phase radiation 
patterns can differ in the way they scatter. 

A. DEFINITION OF RCS 

Radar cross section is a measure of power scattered by the target 
towards the illuminator. The definition of RCS can be stated as 

power reflected to receiver per unit solid angle 
incident power density / 4n 

—► —^ 

In terms of the incident and scattered electric field intensities and Es the 
RCS is 


CT = lim 4nR^ 

R-xx) 



2 


( 1 ) 
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where R is the distance from the target to the observation point. In subsequent 
discussions only the monostatic RCS is considered. That is, the observation 
point is coincident with the transmit source. The limiting process in (1) 
assures that the incident wave is approximately a plane wave. Calculation of 
RCS is essentially a matter of finding the electric field scattered from a target. 
If the current induced on the target by the incident plane wave can be 
determined, then the radiafion integrals used in. noatenna analysis can be 
applied to compute the scattered field. However, determining the induced 
current is a difficult problem because Maxwell’s equations must be solved for 
complicated boundary conditions. 

The unit of RCS most commonly used is decibels relative to a square 
meter (dBsm) which is computed by 

a, dBsm =10 log(a,m^). (2) 

Radar cross section depends on many fectors such as frequency and 
polarization of the incident wave, and the target orientation relative to the 
radar. RCS levels can range from 10,000 m^ for ships to less than 0.001 m^ 
for insects. 

The scattering characteristics of a target are strongly dependent on the 
frequency of the incident wave. There are three frequency regions in which 
the RCS of a target is distinctly different. They are referred to as the low, 
resonance, and high-frequency regions, where low and high are defined 
relative to the size of the target in terms of incident wavelength, rather than 
their physical size. For a target of characteristic length L, the low-frequency 
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or Rayleigh region defined by kL « 1 ( k = In/X). The second is the 
resonance or Mie region where kL » 1. The last is the high-fi-equency region 
or the optical region defined by kL » 1. Depending on the threat radar’s 

wavelength relative to the antenna’s operating band, an arbitrary target can 
fall into either of the three regions. 


B. ANTENNA SCATTERING 

The antenna’s impact on RCS is twofold. First, die antenna fi’equently 
disturbs the continuity of the surface. This introduces more edges with a 
potential for increasing wide-angle scattering. A second aspect of the antenna 
scattering is that the threat radar’s signal can penetrate into the feed and be 
reflected at internal mismatches and junctions. Even though these mismatches 
are small in the operating band, there can be a large number of scattering 
sources that add constructively under some conditions. This scattering is 
dependent on the type of feed and the devices incorporated therein. 

The basic equation of antenna scattering has been presented by several 
authors [Ref 2, 7-9]. The total scattered field for a linearly polarized antenna 
when the port is terminated with a load is 


Es{Zl) = E^iTa) + 


/n 


—> —^ 
h{h 



(3) 


This discussion parallels that given in [Ref 6]. 

6 





where 


Z, = R, + jX, is the radiation impedance 

R, = + Rj the antenna resistance 

Rj = the radiation resistance 

Rj = the ohmic resistance 
—^ 

h = effective height 
—> 

Ei = incident field 

T| = is the intrinsic impedance of the medinm and, 

r 

^ Zl + ZI • 


(4) 


When the load impedance is the complex conjugate of the radiation 

impedance, is called a conjugate-matched load. In equation (3) £t(Z*) is 
the structural mode and the second term on the h^-hand side is the antenna 

mode. When the antenna is conjugate matched, fise antenna mode vanishes 
because the reflection coefficient T goes to zero. This can usually be achieved 
by having real values for the antenna impedance ( i.e., “tuning out” its 
reactance), and connecting the tenninals to a matched transmitter or receiver, 
thereby having no reflection at the junction between the antenna and the 
transmission line. 

Kahn and Kurss [Ref 10] found that for a dipole, conjugate matching 
results in half of the power being scattered and half absorbed. The reason 
being that for an isolated dipole, the structural mode is significant. If low 
RCS is desired, then it is necessary to mismatch the dipole terminals to 
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generate sufficient antenna mode to cancel the structural mode. Forcing the 
antenna mode term to zero by forcing to zero is impractical for large arrays 
because Z, depends on angle and varies with element location in the array. 
However, it is possible in principal to reduce an array’s RCS to low values. 
Green [Ref 11] has shown that if a lossless, linearly polarized antenna 
absorbs more power than it scatters imder matched conditions, its gain in the 
backward direction must exceed the forward direction gain. This condition is 
satisfied in the operating band of high-gain arrays. 

C. FEED REFLECTIONS 

In the case of an operational antenna, the incident threat signal power 
delivered to the radiating element terminals actually enters the feed network. 
The signal transmitted through the aperture elements is reflected at points 
inside of the feed circuit. For a phased array, the feed network can be 
extremely comphcated, with potentially every transmission line discontinuity 
a source of reflected signals. In this thesis, the signal behavior in a series feed 
is examined. 

A series feed is one in which the power to each element is tapped off 
of a main line sequentially, as fllustrated in Figure 1. The coupling values are 
adjusted to provide the desired amplitude distribution. The most common 
method of terminating the main line is matched loading; another is to use a 
short which is called a standing wave feed. For the loaded case, some power 
is always reflected fi’om the load at the end of the main line which gives rise 
to a “reflection lobe” in the radiation pattern. It is reduced in magnitude 
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relative to the main beam by an amount determined by the load reflection 
coeflScient. 

There are two antenna scattering modes to consider in the calculation 
of RCS. The first is the antenna or radiation mode which is determined by the 
radiation properties of the antenna and it vanishes when the anteima is 
conjugate matched to its radiation impedance. The second mode is the 
structural mode that is generated by the induced cmrents on the antenna 
surfaces. This thesis only examines the antenna mode which is the dominant 
RCS contribution to a phased array in the operating band. 

Total RCS can be computed by summing the scattering fi'om all the 
signals that enter the array, are reflected by the mismatches of the 
components, and then return to the aperture. The RCS is obtained fi’om 
combining equations (1) and (3) over all the array elements [Ref 6] 


a(e,<|)) 


4nAe 

~ 




2 


(5) 


where 


r„ (0,(t)) = total reflected signal returned to the aperture for element n 

when the wave is incident fi’om the (0, <}>) direction 

^ = it(ux + ^y + it'z)= wave vector 
w = sin 0 cos ^ 

V = sin 0 sin (|> 
w = cos 0 

d„ = position vector to element n 

= effective area of an element = h^Zg/^Ra 

^norm ~ Dormalized element scattering pattern. 
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There are two approaches to obtaining F in equation (5). One is to use 
a network matrix method such as scattering parameters. All interactions 
between the feed devices are included. The scattering matrix iq)fMroach is 
computationally intense because a matrix equation must be solved. As N 
increases the size of the matrix increases. A second approach is to trace 
signals through the iced However, tracing the signal flow inside a feed 
network can be extremely difiicuh and tedious job, and to include all the 
multiple reflections for accurate calculation is even more difficult. For inband 
cases where the reflection coefficients of the devices are small (r « 1), 
higher order reflections can be neglected The following approximations are 
made for the simplified analysis [Ref. 12]: 

1. Devices of the same type are assumed to have identical electrical 
characteristics. For example, all of the radiating elements have the 
same reflection and transmission coefficients ( and ), with the 
exception of phase shifters which have a transmission coefficient of the 
form 


hn = 


where 


Xn = («-!)«„ 

a, = ^o<f;tsin 05 cos(|», = kod^Us, 

(0,» <l>s) = tfie direction of the array radiation beam. 
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Note that the magnitude of the transmission coefficient is the same for 
all phase shifters ; \tp, | = fp 

2. The feed devices are well matched so that r « 1 and therefore 
hi^er order reflections can be neglected. 

3. Lossless devices will be assumed. 

With the previous assumptions, the fraction of signal entering the 
radiating elements that is reflected back to die aperture for reradiation is 





( 6 ) 


where 


a = kji^u is the incident wave inter-element space delay 
u = sin6cos(t>, 

E'„ represents the signal re-emerging from the series feed at element n. As 
will be seen later, some of this signal is due to elements toward the input (m < 

n\ some due to elements toward the termination (m > n\ and some due to the 
load at the coupler n. 
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Equati(Mi ( 6 ) shows that ^ tc^ severed field is a sum of scattering 
patterns fi^om each mismatch group in the feed 

|£>|=l£; + £J+£J + £;| (7) 

£5 is the total scattered field due to the series feed reflections. It is obtained 
by summing over all as described in the next section. Assuming that only 
one scattering peak dominates at any given incident angle, it may be 

sufficient to approximate the magnitude squared of the sum in (7) by the sum 
of the magnitudes squared 


|£ 1 ^»|£J|^+|£*|%|£J|' + |£5|^ ( 8 ) 

The RCS contribution fi-om the radiating element, phase shifters and 
coupler inputs are [Ref 13] 


_ 47l4^cos^ 6 r sin(JVa) ~|^ 

^ ^ _Nsm{a)_ 


(9) 


Gp 


Gc 


47L4^cos^e 2 fA r sin(Afa) ~|^ 

p *r _^sin(a)_ 


47l4^cos^9 2 fA fA sin(A^) ^ 

Ir ip lNsmiQ_ 


( 10 ) 

( 11 ) 


where ^ = a + a,. 
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D. DERIVATION OF THE SERIES FEED RCS FORMULAS 

Figure 2 shows the signal distribution for the n** four port coupler with 
coupling coefficient c„ and transmission coefficient t.. To begin with, the 



Figure 2: Coupling and transmission paths for a four port coupler. 

couplers are assumed lossless, matched and to have perfect isolation between 
the coupled port and through port. The scattering matrix is 

0 t„ jc„ 0 

fn 0 0 —jCn ^j2'v 

jc„ 0 0 t„ ^ ^ 

0 -jc„ t„ 0 

To obtain formulas for the feed reflections, the received incident threat signal 
can be decomposed into forward and backward traveling waves on the main 
line. Figure 3 shows that a signal received by element n appears as an 
excitation for the N - n elements toward the load and the n - 1 elements 
toward the input. Other scattering terms not included in the forward or 
backward beam are the “self-reflected beam” due to the reflection from the 
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load terminating element n, and the *'input load reflected beam” due to the 
reflection from the input load or receiver at the beginning of the feed. Thus 
for each element in the array four beams are reradiated, and their beamwidths 
and peak levels depend on the element location in the array. The total 
scattered field is obtained by the superposition of the forward, backward, 
self-reflected, and input load reflected beams from all elements. 


input port 
loaded 

< -► 






Figure 3 : Forward and backward traveling waves in a series feed network. 

The formulas for the RCS of the series network can be obtained by 
tracing signals through the feed. Let be the reflection coefficient of all of 
the loads. For an incident signal at element n, “forward wave” that travels 
toward element N (see Figure 4), radiates reflected signals from n+1 up to 
element N and can be expressed as 


(2) Tbe expressiOQs for the electric field intensity (E) have the leading constants and free space 
Green's function omitted because only a ratio is needed in Equation (1). 
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Ef^ = eA"-lXa,4a) { e>'*‘jCn*i e'"<“'-~) 

+ ^'0- t„HjCn*2 + ... 


+ /«l/«2 • - tN-l JCN ^xr-ixa,™) } _ (13) 

where Vo~kL. The electrical line length between couplers is denoted as 
and it is equal to k times the physical length of the line. The “backward 
wave” that travels toward the input (see Figure 5), radiates reflected signals 
from the n-1 elements downline from element n 



Figure 4 : Forward wave (traveling toward element N). 


+jc„-\e ^'^-^Ti ° (oa-Kx)j 
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The self-reflection is due to the load at coupler n, and it is unique in its form. 
Thus it can be separated and simplified as 



Figure 5 : Backward wave (traveling toward element 1 - input). 

For the reflection from the input load, the equation is given by 

n ( 16 ) 

The total forward, self, input load, and backward traveling fields from 
all the elements caused by an incident wave is obtained by summing over all 
elements 
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^-1 ( 

s 

f m-1 

(17) 

M-l 1 


V|=1 J 

d 

N-\ 

r n-l 

fn-\ \ ] 


Eb = -TL Z 

I C 


(18) 

h^2 

L»»»i 

V|«m / J 



(19) 

E„i/^ f Ti tl 

■ 


(20) 


The minus sign occurs for some terms because the reflected beams pass 
through two coupled arms yielding a f factor. The total scattered field due to 
the series feed beam is 

|£ 5 |' = \Ef + Et+ E,„ + 

’«\Ej\^ + \Ei\^ + \Ej^ + \Esei/ ( 21 ) 

The correct normalization can be obtained by comparing the total 
scattered field in equation (21) to the maximum that would be obtained by an 
array with completely reflecting elements (rN). Thus 





The coupling coefBcients (cj are determined by the prescribed radiation 
beam sidelobe level. Appendix A describes how the coefficients are 
computed from a given amplitude distribution. 
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III. COMPUTATION AND RESULTS 


Id this chapter, calculated results are presented for series fed arrays 
using equations (17) through (22). To verify the approximate method, data is 
compared with the rigorous method using scattering matrices. The scattering 
matrix solution is described briefly in Appendix C. The aperture element is 
represented by a two-port device with the reflection coefficient allowing 
direct comparison between both approximate and scattering matrix solutions. 
This provides a means of evaluating the validity of the assumptions made in 
deriving the approximate equations. 

For the series-fed phased array with series feed netwcHk operatmg at 
wavelength X, calculations were performed for the parametn^ sets listed in 
Table 1. 


TABLE 1 : SERIES-FED PHASED ARRAY PARAMETERS. 


parameter 

range 

no. of elements, N 

50 

Vo 

0, 7 t/ 4 , 7i/2, n 

d 

0.4 X 


0.2 

1 

0.5 X 


uniform, 

cosine-squared on a pedestal 

e. 

0^45‘’ 
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A MATLAB program for calculating the approximate RCS is shown in 
i^pendix B. Each element of the array (Figure 1) that is illuminated by the 
incident wave has its own forward wave, backward, self and input reflected 
term as discussed in Chapter H, Section D. 

Some parameters are varied to study the behavior of RCS. They 
include the scan angle (0^), the line length between couplers ( ), and the 

coupling coefficients ( c„). Figures 6 and 7 show the difference between the 
broadside and scanned RCS when \|/^ = 7c/4. It can be observed that the 
highest spike at 6 = 0“ on both grajrfis is due to specular scattering from the 
radiating elements. The spikes at about >20® in Figure 6 and at 20® in Figure 7 
are due to the reflection from the input load. The spike at -35® in Figure 7 is 
due to Bragg diffraction, and the one at 45® is due to the effects of forward 
and backward waves on the main line scattering retrodirectively. A 
comparison of the scanned and broadside cases shows that the reflections 



Figure 6 : RCS of series fed array by the approximate method. 
(0s = O®, v. = «/4) 
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from points behind the phase shifters move toward or away from the specular 
beam along with the array’s radiation beam. 



Figure 7 : RCS of series fed array by the approximate method. 
(0s = Jt/4, = tc/4) 


In general, the individual spikes can be attributed to specific scattering 
sources in the antenna. The specular reflection ( due to r,) has its peak at 9 = 
0°. There are no Bragg lobes for this term because the element spacing is 
sufficiently small. Reflections from the input of the phase shifter have a 
similar behavior as can be seen from equation (10). Reflected signals from the 
coupler inputs will pass through the phase shifter and therefore scan along 
with the radiation beam. (This can be a problem for an array on a stealthy 
platform if it is used to track a threat.) When the array is scanned, a second 
lobe due to Bragg diffraction moves into the visible region. 
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Figure 8 : RCS of series fed array by the approximate method. 
(0s = 0®, V, =n/2) 



Figure 9 : RCS of series fed array by the approximate method. 
(0s = it/4, V, = ii/2) 
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A large lobe also arises due to the reflection from the input load. In this 
case, the array acts as a receive antenna that collects the incident threat 
signal. A portion of this signal is reflected at the terminals and appears just as 
if a transmit signal has been injected into the antenna. Thus the RCS pattern 
is the gain pattern squared. It can be seen by comparing Figures 6, 7 and 9 
that the lobe from the input load reflection is determined by \\i^ as well as 0,. 



Figure 10 : RCS of a series fed array by the scattering matrix method. 

(0. = O*, V. = Jt/4) 

The RCS patterns obtained using the scattering matrix method are 
shown in Figures 10 through 13. They correspond to the same cases shown in 
Figures 6 through 9. This solution can be considered rigorous in the sense that 
all interactions and multiple reflections are included. The agreement between 
the approximate and the rigorous solution is excellent. 
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Tbeta(deg) 


Figure 11 : RCS of a series fed array by the scattering matrix method, 
(0, = 7C/4, V. = 5C/4) 



Figure 12 : RCS of a series fed array by the scattering matrix method. 

























Figure 14 : RCS of a series fed array with tapered amplitude distributiou by 
the approximate method (6,«ic/4, y, = n/2 ). 
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Figure 15 : RCS of a series fed array with tapered amplitude distribution by 
the approximate method ( 6, = 0*, \|f, = ji/2 ). 

The effects of changing the coupling distribution is illustrated in 
Figures 14 and 15. The RCS behavior is essentially the same but the feed 
contributions to RCS are reduced slightly. 
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IV. CONCLUSION 


i^proximate formulas for the inband RCS of an array with series feed 
have been derived. The formulas are based on the hypothesis that an incident 
wave excites forward and backward traveling waves on the main line. The 
approximate formulas are in good agreement with results obtained using 
scattering matrices, thereby verifying die assumptions made in the 
approximate solution. 

Spikes in the RCS pattern have been identified with specific scattering 
sources in the array. The parameters afifectiiig the level and location of the 
lobes have been noted. By properly choosing parameters such as \j/„ and T^, it 
is possible to control the spikes to some extent. The data shows that there is 
no preferred value of \\f^ that eliminates the input load reflection. It is 
common practice that the line length between couplers should be set to an 
odd multiple of a quarter wavelength to cancel multiple reflections. However, 
total suppression of the RCS can only be achieved by perfectly matching all 
of the feed devices, which is practically impossible. 

The main advantage of the approximate method relative to the 
scattering matrix method is its speed. The scattering matrix method requires 
that a matrix equation be solved, and its size increases with the number of 
array elements. A typical approximate calculation takes about 2 mm. on a 
Sparc 10, whereas the matrix solution runs about 2 hours (50 elements and 
0.5° increments ). Also, the approximate formulas are easily extended to an 
arbitrary number of elements. Since both methods assume lossless feed 
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networks and identical reflection coefScients for all similar devices, the actual 
value of RCS can be expected to be slightly lower than shown. 

For future research, it is recommended that the RCS of a 
two-dimensional array be investigated and be presented in contour form [Ref 
6]. The RCS of hybrid feed networks (e.g., series fed in one dimension and 
parallel fed in the other) should also be examined. 
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APPENDIX A : COUPLING COEFFICIENTS 


This appendix derives the series feed coupling coefficients for a given 
amphtude distribution {o„}. For uniform illumination (13dB sidelobe level) 
=1, n=l,2, N where N is the number of couplers. Also of interest is a 
cosine-squared on a pedestal distribution (32dB sidelobe level) where 



Figure A.l : Cosine-squared distribution. 
an =ao +(I- Qo) sin^(^ x 7t) 

and Go = 0.2, as illustrated in Figure A.I. Referring to Figure A.2, for a 
lossless feed 

Pin =11 kal+Pi 

n=l 

or. 
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Pin-Pi^k'Lal 

n 


and 


r. _ f«-n fgr ‘ l-Pl 
v ,.2 ^ v ^2 


ILo 

n 


Note that is the power into the load and is the total input power. 



2 2 

kfl, 




-VV\A 


Figure A.2 : Power into the lossless feed. 


The coupling coeflScients are defined by 

"" p 

m 


where 


Pm 



m-l _ m-1 _ 

2 kal = 1-2 kal 

q=\ ^ 9=1 ^ 
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Therefore, for unit power into the load (P^ =1), the coupling coefficient 
becomes 





m-\ 
Pin- Z 
q=\ 



N «-l 

Tk 2 “i- s 


n=l w=l 


The transmission coefficients are given by 
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APPENDIX B : MATLAB PROGRAM LISTING 


This appendix presents a listing of the MATLAB Code that computes 
the RCS using the approximate formulas. The version hsted is for uniform 
amphtude distribution and no scan. 

%% INBAND SCATTERING FROM ARRAYS WITH SERIES FEED NETWORKS 
%% Approximate Method (Broadside RCS Uniform Distribution). 

N=50;, Gamma=0.2;, Si=pi/4; 

Ao=l;, Pl=0.05;, ki=0.8*pi*sm(0); 

theta=-pi/2:pi/360:pi/2;, thetal=theta.*l 80/pi; 

% a) Forward Wave ; (Travels toward element N) 

for a=l;N; 

A(a)=Ao+((l-Ao)*(sin(((a-l)/(N-l)).*pi)).^2); 

end 

for 1=1 :N; 

Cm(l)=sqrt((A(max(l))^2)/(((l/(l-Pl))*(sum(A(l:N).^2)))-(sum(A(l:l-l).^2))) 

); 

end 
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for i=l;N-l; 

tl (i)=((sqrt(l -(Cm(max(i)).''2))) *exp(j*Si)); 
t2(iHsqrt(l -(Cin(max(i)).^2))); 


fork=l:N; 

tsf(k)=(l -(Cm(max(i)).^2)); 
t3(k)=prod(tl(l:k-l)); 
t4(k)=prod(tl (1 :k-2)); 
t5(k)=prod(t2(l ;k-l)); 


for z=l:361; 

del 1 (z)=sin(theta(max(z))).*(0.8*pi); 

del=dell(max(z)); 

for n=l;N-l; 

for m=n+l:N; 

Ef5(m)=Cm(m).*exp(j*(m-l)*(ki+del))*t3(m); 

end 

Ef4=sum(Ef5); 

Efi(n)=Cm(n). *exp(j*(n-1 )*(ki+deI))*Ef4; 
end 

Ef2=sum(Ef3); 

Efl (z)=-Gamma*Ef2; 
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% b) Self term ; 


forp=l:N; 

Esf3(p)=Gamma*tsf(p) ♦exp(j*2*(p-l)*(ki+del)); 
end 

Esf2=sum(Esfi); 

Esfl(z)=(Es£2); 

% c) Backward Wave : (Travels toward element 1-input) 


for q=2:N-l; 

for r=l:q-l; 

Eb5(r)=Cm(r).*exp(j*(r-l)*(ki+del))*t4(r); 

end 

Eb4=sum(Eb5); 

Eb3(q)=Cm(q).*exp(j*(q'l)*(ki+del))*Eb4; 

end 

Eb2=sum(Eb3); 

Eb 1 (z)=-Gamma*Eb2; 

% d) Input load reflected wave 

for s=l:N 

Em3(s)=Cm(s) ♦exp(j*(s-l)*(ki+del+Si))*t5(s); 
end 
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Ein2=sum(Ein3); 

Einl (z)=-Gamma*(l *Gamma^2)^2*(Em2).'^2; 
end 

Ef(:)=10.*logl0(abs(Efl(:))); 

Estl[:)=l O.*logl 0(abs(Esn (:))); 

Eb(:)=10.*logl0(abs(Ebl(:))); 

Em(;)=l O.*logl 0(abs(Eml (;);); 

Es=(abs(Efl)).^2+(abs(Esfl))."'2+(abs(Ebl)).^2+(abs(Einl)).'^2; 

Enl=Es./(N'^2); 

for s=l:361; 

if sin(dell(max(s)))—0 
Rf(s)=l; 
else 

Rf(s)=(sin(N*del 1 (max(s)))./(N.*sm(dell (max(s))))).^2; 
end 

if sin(0.8*pi*(sin(theta(max(s)))+siD<'0)))=0 
Rfs(s>=l; 
else 

Rfs(s)=(sin(N*0.8*pi*(sin(theta(max(s)))+sin(0)))./(N. ♦sin(0.8*pi*(sm(theta( 
max(s)))+sm(0))))).^2; 

end 

end 
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fort=l:361; 

RCSl(tH400*pi*(cos(theta(max(t))))^2).*(l-Gamma^2r4.*Enl(t); 

RCSr(t)=(400*pi*(cos(theta(max(t))))^2).*Gamma'^2.*Rf(max(t)); 

RCSp(t)=(400*pi'''(cos(theta(max(t))))'^2).*Gamma^2*(l-Gamma'^2)'^2.*Rf( 

max(t)); 

RCSc(t)=(400*pi*(cos(theta(max(t))))'^2).*Gamma^2*((l-Gamma'^2)^4.*Rfs 

(max(t))); 

end 

RCS2(;)=10.*logl 0((RCS 1 (:))+1 .e-10); 
RCSR(;)=10.*logl0((RCSr(;))+l.e-10); 
RCSP(:)=10.*logl0((RCSp(;))+l.e-10); 
RCSC(:)=10.*logl0((RCSc(:))+l.e-10); 

RCSt=RCS 1 +RCSr+RCSp+RCSc; 

RCS(:)=10.*logl0((RCSt(:))+l.e-10); 

figure 

plot(thetal,Ef), grid, axis([-90 90-50 30]); 
title(Torward wave vs Angle') 
xlabel('Theta(de,g.)'), ylabel(Ef(dB)') 
pause 


36 










figure 

plot(thetal,Esf), grid, axis([-90 90 -50 30]); 
title('Self-Reflected wave vs Angle') 
xlabel(Theta(deg.y), ylabel('Esf(dB)’) 
pause 

figure 

plot(thetal ,Eb), grid, axis([-90 90 -50 30]); 
titleCBackward wave vs Angle') 
xlabel('Theta(deg.)'), ylabel('Eb(dB)') 
pause 

figure 

plot(thetal,Ein), grid, axis([-90 90-50 30]); 
titleCLoad reflected wave vs Angle') 
xlabel('Theta(deg.)'), ylabel(Ein(dB)') 
pause 

figure 

subplot(221), pIot(thetal ,RCS2) 

grid, axis([-90 90 -50 30]);, title('Ant. RCS vs Angle*) 

xlabel('Theta(deg.)'). ylabel(TlCS(dB)') 

pause 
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subplot(222), pIot(thetalJlCSR) 

grid axis([-90 90 -50 30]);, titleCReflected RCS vs Angle') 

xlabel('Theta(deg.)'), ylabel(RCS(dB)’) 

pause 


subplot(223), plot(thetal,RCSP) 

grid axis([-90 90 -50 30]);, title(Phase Shifter RCS vs Angle*) 

xlabel('Theta(deg.)'), ylabel(RCS(dBy) 

pause 


subplot(224), plot(thetal,RCSC) 

grid axis([-90 90 -50 30]);, title('Coupler RCS vs Angle') 

xlabel('Theta(deg.)'), ylabelCRCSCdB)') 

pause 


figure 

plot(theta 1 ,RCS), grid 

axis([-90 90 -50 30]);, title('Total RCS vs Angle’) 

xlabel('Theta(deg.)'), ylabel('RCS(dB)') 

save rcsbn2 thetal RCS Ef Esf Eb Ein RCS2 RCSR RCSP RCSC 
pause, axis;, end 
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APPENDIX C : SCATTERING MATRIX FORMULATION 


The series feed is represented by an interconnection of two and four 
port devices as shown in Figure C.l. The unknown quantities are the set of 
A scattering equation can be written at each port. Thus for an N 
element array, there are %N equations that must be solved simultaneously. 

The scattering matrix of the radiating element is 


' dn dn ■' 


rr tr 

_ d2\ d22 _ 


_ir rr _ 


For the phase shifter 

rp ‘ 

The excitations are due to the incident plane where S, = . 

Starting with element 1, the two scattering equatioos at the radiating devices 
are 


Pn = 


(P")n ^>^12 
I (Pn)2i iPn)22 


a\=S\d\i-^a3d\2 ( 1 ) 

02 = S\d2\ + ( 2 ) 
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and at the phase shifter ports 


as = a2(pi)ii +^5(Pi)i2 (3) 

04 = a2(pi)2i '^^^(P022 

Four equauoiis can be written at coupler 1 (cj 

06 = 00 r(ci)ii +V|/af 9 (ci)i 2 +i 3 f 4 (ci)i 3 +flr 7 r(ci)i4 ( 5 ) 

Oi = 06 r(ci) 2 i +M/ a9(ci)22 + 04 ( 01)23 + a? r(ci)24 ( 6 ) 

<35 = 06 r(ci)3i +11/09(01)32+04(01)33 + 07 r(ci)34 ( 7 ) 

07 =06 r(Ci)4i +V|/flf9(Ol)42+04(01)43+07 ^(01)44 (8) 

It has been assumed that all loads are identical (r= r,= ... = F). 

Similarly, all interelement line lengths are equal (h/i= ~ 

The coupling coefficients are given in Appendix A and the scattering matrix 

by equation (12). 

For elements 2 through N, a pattern repeats. At the radiating element 

08n-6 = Sn <^11+ 08«-4 d\2 (8(n-l) + 1) 

08n-5 = Sn 0^21+ 08n-4 d22 (8(n-l) + 2) 

and at the phase shifter 

08«-4 = 08n-5(p«)ii + 08n-2(Pn)i2 (8(n-l) + 3) 

08m-3 = Oi„-5(pn)2i + 0%n-2(Pi^21 (8(n-l) + 4) 
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For coupler 2 through N-1 


a%n-l = w ai„~i(Cn)ii + V|/ a%n+liCn)i2 + ^8«-3(^^«)l3 

+ai„-\r(cn)i 4 (8(n-l) + 5) 

ain = Vj/ ain-sic„) 2 i + \\f as„+l{c „)22 + fl8»-3(^^«)23 

+ain-iTic „)24 (8(n-l) + 6) 

ain-2 = M/ ain-i(Cn)3\ + M/ a^n+l(Cn)32 + ^Sn-siCrds^ 

+ai„-\T(c „)24 (8(n-l) + 7) 

Oin-l — ^8«-8(^«)41 + OirH-l(Cn)42 + ^in-3(Cn)43 

+<38«-ir(c„)44 (8(n-l) + 8) 

Finally, coupler N must be treated separately because of the extra load 


a%N-7 = M/ ^8A^-8(^A^)ll + ^SN T(cn)i2 + <^iN-3(CN)i3 

+a iN -1 r (cn) ] 4 (8N-3) 

azN = V ar87V-8(CA^)21+ r(cA?)22 + 0%N-3(Cn)23 

+ci%ff-\T(cN)24 (8N-2) 

OiN-l = v aiN-i(cN)31+ Cl%N ^(Cn) 32 + a^N-3{CN)33 

+CI%N-\^{cn) 34 (8N-1) 

= V ^3r8A^-8(CAr)41+ Cl%N r(CAr)42 + a%N-3{CN)43 
+CliN-\^{CN)44 (8N) 


The above 8N equations can be cast into matrix form 


[S\ = [F\[a] 






where the vector [S] contains the excitations {S„} and [a] the unknowns {a„}. 
[F] is a feed scattering matrix obtained from the coefficients of the above 
equations. They can be represented as entries in the following table 


column = unknown index 



The matrix equation is solved to obtain and the reflection coefficients r„ 
in equation (5) can be calculated from ai, ajo,flgw-e, •••, ci%n-6- 
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Figun C.1 : Scattering matrix feed diagram. 
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